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Recent observations from PAMELA, FERMI, and ATIC point to a new source of high energy 
cosmic rays. If these signals are due to annihilating dark matter, the annihilation cross section in 
the present day must be substantially larger than that necessary for thermal freeze-out in the early 
universe. A new force, mediated by a particle of mass 0(100 MeV), leading to a velocity dependent 
annihilation cross section - a 'Sommerfeld enhancement' - has been proposed as a possible expla- 
nation. We point out that such models necessarily increase the dark matter (DM) self-scattering 
cross section, and use observational bounds on the amount of DM-DM scattering allowed in various 
astrophysical systems to place constraints on the mass and couplings of the light mediator. 

INTRODUCTION 

It is well known that the Standard Model (SM) particles are not the primary constituents of the Universe. Dark 
energy makes up 72% of the energy budget, while dark matter makes up 24%, and only the remaining small fraction 
is due to the SM [T] . The nature of dark matter (DM) is one of the major puzzles facing physics today: other than 
its energy density and gravitational interactions, little about it is known with any certainty. 

Many solutions to the hierarchy problem naturally contain new stable massive particles. Their masses are typically 
at the weak scale and their annihilation cross section is of the right order of magnitude such that, if they were thermally 
produced in the early universe, they would have the correct abundance today to be the DM. For this so-called "WIMP 
miracle" (Weakly Interacting Massive Particle) to occur, and the DM to be a thermal relic, the annihilation cross 
section at the time of freeze-out, when the temperature was T ~ m x /25, must have been (<jv) ~3x 10~ 26 cm 3 /s. It 
is usually assumed that this annihilation cross section is velocity independent and that in the present epoch the DM 
is annihilating with the same cross section. This is an appealing possibility, but by no means the only one. 

Recently there have been a series of observations of high energy electrons, positrons, and gamma rays from HEAT 
0, AMS-01 0, PPB-BETS gj, PAMELA [5], ATIC 0, and the Fermi Gamma Ray Space Telescope (FGST) 0. 
These measurements are seemingly at odds with what is predicted from the secondary production of anti-particles 
from cosmic-ray propagation. These excesses may be explained by astrophysical processes - for instance nearby 
pulsars [SI [S] or supernovae remnants |10j may be a source for high energy positrons and electrons - or they may be 
due to DM annihilating in our galaxy. 

If these anomalies are related to DM annihilation 11-26], the necessary value of (av) in the present day must be 
larger than the thermal freeze-out value of (av) ~3x 1CU 26 cm 3 /s by a factor of approximately 10 2-3 (see e.g. 27]). 
It is possible that this enhancement is purely structural in nature: if sufficient halo substructure exists wherein the 
DM density, p, is larger than that expected from Galactic simulations, the annihilation rate, which scales as p 2 , would 
be increased. Alternatively, the presence of a new force interacting with DM and mediated by a light boson (scalar 
or vector) could provide a velocity-dependent annihilation cross section via the Sommerfeld effect [TTJ [2"5H5U] . This 
enhancement typically scales as i> -1 . However, due to these new forces the DM particles x can develop bound states, 
near which the enhancement is much larger, proportional to ~ v~ 2 . 

For the remainder of this paper, we shall take the Sommerfeld enhancement via some massive (but relatively light) 
boson <j> as the explanation for the cosmic ray data. We then are interested in the associated enhancement of the 
self-scattering cross section of dark matter (i.e. \~X or X~X scattering, not x-nucleon scattering). The same non- 
perturbative effect that alters the wavefunction at zero distance will also cause a phase shift at infinity. As a result, 
an increase in the scattering cross section goes hand in hand with the Sommerfeld boost to the annihilation rate. 
However, as we shall see, though both these boosts occur at the same points in parameter space, the numerical size 
of the boost differs greatly between scattering and annihilation. 

It has been shown, both by A^-body simulation and direct observation of the system of colliding galaxy clusters 
known as the Bullet cluster, that dark matter must be, to good approximation, collisionless. From these observations, 
we may place limits on the mediator mass and coupling of any new force at work in the dark sector. Constraints from 
dark matter collisions in the case of a massless mediator have been considered previously [311 [22] . 

In the remainder of this paper we first derive a general expression for the scattering cross section of two dark 
matter particles interacting via some light force carrier. This expression must be solved numerically, so we also find 
an approximate formula for the cross section, which is valid in the regime of very small mediator mass. After this, 
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we consider the various bounds on the self-scattering cross section and the characteristics of the systems from which 
these bounds are extracted. Using the approximate parameters of these systems we then constrain the coupling and 
mass of possible new dark forces. 



SCATTERING AND ANNIHILATION 



We wish to consider the interactions between two slowly moving DM particles exchanging light force carriers. As 
the DM particles are non-relativistic they may exchange multiple bosons while undergoing either an annihilation 
or scattering process, as shown in Fig. [I] This is a non-perturbative effect that must be resummed, which is done 
by solving the Schrodinger equation for the reduced system. The potential between the two particles alters the 
wavefunction of the reduced system both at r — 0, affecting the annihilation rate, and at r — > oo, affecting the 
scattering cross section. Both scattering and annihilation involve similar diagrams, the only difference between the 
two is an insertion of the short distance operator responsible for DM annihilation. Thus, whenever the parameters of 
the system are such that there is a large enhancement in the annihilation cross section, then the self-scattering cross 
section is also enhanced. Since the annihilation diagram must end in the x~ X self-annihilation interaction which is 
not present in the scattering, there is no reason to expect that the size of both enhancements will be equivalent. 

We consider two DM particles of mass m x , with the force between them given by a potential V(r). For the rest of 
this article we will restrict ourselves to the case of a Yukawa potential, generated by the exchange of a boson of mass 
which couples to the DM with coupling strength A: 



V(r) = ~- e - m * r . 
r 

Here a — A 2 /An. The wavefunction of the reduced system is i/j(r) 
function, Rp(r), satisfies the radial Schrodinger equation, 



(1) 



E<m^( r )^,m(^)' where the radial wave- 



1 d 

r 2 dr 



j dR( 
dr 



k 2 



£(£+!) 



2/jV{r) 



Ri = 



(2) 



Here, /U = m x /2 is the reduced mass and k = /iv Te \ is the momentum in the reduced system. It is useful to introduce 
Xe = rRt and x = m x ar, in terms of which Eq. ^ becomes 
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We solve this equation with the boundary conditions that x is regular at the origin and at large r behaves as 
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The differential scattering cross section is given by 
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and the total scattering cross section is then given by the sum over all angular momenta I: 
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A quantity of interest when discussing observational constraints is the transfer cross section, cr t r, which is a weighted 
average of the differential cross section 
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FIG. 1: Non-perturbative 
scattering (right). 



"ladder" diagrams corresponding to the formation of a \ ~ X bound state for annihilation (left) and 



This weighted cross section controls the rate at which energy is transferred between particles in a collision. 

For the case of annihilation the enhancement is determined from the value of the wavefunction at the origin, 
|^(0) | 2 . At low velocities the attractive potential distorts the wavefunction, increasing it at the origin - the Sommerfeld 
enhancement [H]. For a Yukawa potential this enhancement at low velocities scales as ~ 1/v 29, 33, 34J but saturates 



at velocities of order v ~ y 'am ( p/m x . For some low velocities there are particular points in parameter space where 



there is a light resonance due to a bound state in the potential [35H37] . which can greatly increase the annihilation 
cross section. These same bound states will also lead to a large enhancement in the scattering cross-section. 

As the Yukawa potential does not have analytical solutions, we are required to numerically solve Eq. Hence, 
we cannot sum over all possible values of the orbital angular momentum £ to find the total cross section. It is 
therefore useful to determine what the largest value of £ that is relevant to the sum in Eqs. (|6| and ([7|. We can 
estimate the relevant angular momentum of the system by considering the momentum k, and impact parameter b, 
of the incoming particle. We would expect L, the largest angular momentum needed to describe the interaction, to 
be L = fc& max = fJ,v Te ib max . We estimate & maX j the largest impact parameter relevant to the scattering, to be the 
separation at which the potential energy is comparable to the kinetic energy. The potential is finite range, so for 
t ^S> &max there is little scattering and the contribution of the corresponding £-modes to the total cross section should 
be small. 

For values of large relative to the DM's kinetic energy, it is clear that 6 ma x 5s ■> since the potential is 
controlled by the exponential. In general the largest impact relevant impact parameter is given by the solution to the 
equation, 



1 2 a 
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The resulting values of L are shown in Fig. [2] for two values of the coupling: a = 0.1 and 0.01. 

Before solving for the phase shifts numerically it is useful to consider a simple approximation of the cross section in 
a Yukawa potential. We make the approximation that for all £ < L the phase shift is maximal i.e. sin Sg ~ 1 and that 
for all £ > L, the phase shift is zero. This approximation is borne out in numerical simulation, when L > 1. Under 
these assumptions, the total cross section, Eq. Q, can be approximated as 
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If the phase shifts are exactly maximal, sin<5f = 1, then there is considerable cancellation that occurs in the sum for 
the transfer cross section In this limit it is 
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However, there is unlikely to be exact cancellation between the different phase shifts so one would expect, for the 
case with multiple £- modes contributing, that the answer lies between ^ and (10). Note that in the limit that the 
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FIG. 2: The maximum angular momentum L = m x u6 max as a function of the mediator mass for DM of velocity v colliding 
head-on. In the left-hand plot we have chosen a = 0.01 and in the right-hand plot a = 0.1. In both cases m x = 500 GeV. Each 
contour line corresponds to an increase of L by one from the previous. 



potential is turned off, a — > 0, V{r) can be treated as a perturbation and the Born approximation is reasonable. Our 
approximation is not valid in this limit as we assume that the scattering is maximal for all I < L, whereas in the Born 
limit L = and i5o ~ am x k/rr&, making the cross section velocity independent. 

We explicitly calculate the transfer cross section by summing the phase shifts for the low-lying ^-modes. The region 
of parameter space for which there are observational constraints has velocities in the range 10 km/s < v < 1000 km/s, 
and we limit ourselves to mediator masses above ~ 10 MeV. In Fig. [3] we present the resulting <7t r , for DM of fixed 
speed of 100 km/s and two choices of a (a = 0.01 and a = 0.1). For < 500 MeV we sum the contributions from 
I < 5 and for larger only the first two modes are included. In Fig. [4] we illustrate the dependence of 0t r upon v for 
fixed mediator mass, for the same two values of a. Although we have summed up several i modes, for all but very low 
mediator masses the cross section is dominated by the s-wave. As anticipated, there are points in parameter space 
with nearby bound states, causing a resonance in the scattering which also results in an increase in the annihilation 
rate (dashed blue line in Fig. |3| . 

It can also be seen in Fig. that our simple approximation (dotted red lines) does a good job of capturing the 
correct behaviour at low mediator mass, but at higher mediator masses it significantly overestimates the result away 
from the positions in parameter space where there is a resonance. This is as expected, since at a resonance sin 5# — > 1, 
which is the saturation assumption made in the derivation of the approximate formula Eq. (10). 

This approximation is also useful since solving the Schrodinger equation ^ numerically for a large number of 
i'-modes is time consuming. As mentioned earlier, if the coupling, a is weak or the velocity of the DM large then 
the potential is a small effect and the Schrodinger equation can be solved perturbatively, the Born approximation. 
However, for the low relative velocities of interest here the Born approximation is not valid over much of the parameter 
space. Alternatively, one can numerically solve the classical motion of a particle in the same Yukawa potential. For 
a large range of parameters these numerical results are well fit by a simple analytic expression [38| 139] . While the 
classical result will miss the existence of resonances present in the quantum mechanical solution, since classically there 
can be no tunneling, the two approaches agree away from these special, but not negligible, regions of parameter space. 



OBSERVATIONAL BOUNDS 



Self-interactions transfer energy between DM particles and thus tend to drive DM halos into a spherical isothermal 
configuration, with an approximately Maxwellian phase space density. The observation of non-spherical halos and 
halos with large phase-space densities (dwarf spheroidals and the cores of galaxy clusters) place constraints on the size 
of the self-interactions. A weaker bound has also been derived from the Bullet Cluster. These bounds are summarized 
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FIG. 3: For m x = 500 GeV, and a — 0.01 (left-hand plot) and a — 0.1 (right-hand plot) we show the transfer cross section, 
the numerical results (blue solid line) and our approximate formulae Eqs. |9| and (10 1 (upper and lower red dotted lines), as 
well as the Sommerfeld enhancement (blue dashed line) in the annihilation cross section. We have assumed that the DM is 
colliding head-on with speed 100 km/s. 
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FIG. 4: For m x = 500 GeV, — 100 MeV and a = 0.01 (left-hand plot) and a = 0.1 (right-hand plot) we plot the transfer 
cross section for two DM particle colliding head on at speed v. The result of the numerical calculation, summing the first five 
^-modes, is shown in blue (solid) and the upper and lower red (dotted) curves uses the approximate cross section Eqs. |9| and 
( 10 1 described in the text. 



in Table I, along with the characteristic velocities of the relevant systems. We review them individually here. 1 

The Bullet cluster bounds [4D - I43] are weak, but conceptually simple. The Bullet cluster system consists of two 
galaxy clusters which have collided. The baryonic gas has been slowed by the collision, while gravitational lensing 
demonstrates that the two DM halos have flowed freely through each other. From this, a robust bound of cr/m x < 
1.25 cm 2 /g has been placed on the self-scattering; however, by making additional assumptions about the initial states 
of the two clusters one can strengthen this bound to < 0.7 cm 2 /g |43) . 

Smaller DM halos must be colder than larger ones, as fast moving DM particles can escape from a smaller grav- 
itational well. If a smaller subhalo is in orbit around a larger halo (i.e. dwarf galaxies around a spiral galaxy or a 
spiral galaxies around a galactic cluster) then large cross sections allow the efficient transfer of heat from the hot DM 
in the large halo to the colder DM in the subhalo. This causes the subhalo to dissipate, either through spallation or 
evaporation (depending on the size of the scattering cross section) [33] . In order that galaxies such as our own do not 
evaporate too quickly in the hot DM of the host cluster, the range 0.3 cm 2 /g < a/m x < 10 4 cm 2 /g may be excluded 
[45] . In addition, if the cross section is allowed to depend on powers of velocity, i.e. cr/m x = (To/rn x (v/v )~ 2S , then 



1 The unit of choice for <r/m x in TV-body simulations is the cm 2 /g. We adopt this convention for the remainder of the paper. For 
convenience, note that 1 GeV -3 ~2x 10 -4 cm 2 /g. For an overview of the constraints, see [40] . 
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the rate of evaporation of dwarf galaxies rules out the range 1.1 < 5 < 1.8 for aa/m x = 0.1 cm 2 /g [35J. For such 
bounds the characteristic velocity is that of the larger, hotter halo, rather than the smaller and colder one. 

In addition to causing the dissipation of smaller halos, self-scattering also causes the cores of DM halos to become 
more spherical, larger, and less cuspy, as the increased rate of energy transfer allows smoothing over the central 
region. As the cores of some galaxy clusters are known to be elliptical f4"6l!48| . a limit of cr/m x < 0.02 cm 2 /g can be 
derived from these systems [35]. Previously, the lack of cusps in the cores of dwarf galaxies, combined with a paucity 
of dwarf subhalos relative to the number of large galaxies in the Local Group (the 'missing satellite problem'), was 
taken as evidence in favor of collisional dark matter [49]. The range of 0.5 cm 2 /g < cr/m x < 5 cm 2 /g was claimed 
to put TV-body simulations in better correspondence with observation [331 |4"5H53"] . However, it since has been shown 
that tidal stripping |54j or photoionization effects [55) can bring simulation in line with observation. In this case, 
no additional scattering is necessary to explain the structure of dwarf galaxies, and indeed, we can place a bound 
of cr/m x < 0.06 cm 2 /g by requiring that the TV-body predications do not differ from collisionless DM [55]. To be 
conservative, we take this dwarf galaxy bound to be 0.1 cm 2 /g. 

Studies of the thermodynamics of galaxies indicate that if the average time between collisions is less than Hubble 
time, then the cusps of dark matter cores would be flatter and larger than what we observe [53] [57J [55] ■ In order that 
this not occur, the cross section is be bounded by [57] 

a . , . / 0.02 Af pc 3 \ / 100 km/s\ 

< 0.2 cm 2 /g ^— 1- . (11) 

m x \ p J \ v Q J 

Here, p is the DM density of the system. The DM systems considered in Ref. [57] have velocities and densities such 
the cross section is limited to be < 0.01 — 0.6 cm 2 /g. 

Finally, the rate of growth of super-massive black holes (such as the one in the core of our own Galaxy) place a 
limit on the scattering of dark matter [59j [60] . If dark matter is collisional, then supermassive black holes would 
grow faster than in the collisionless case; rapidly reaching the point where the size of the accretion disk approaches 
thcmean free path of the DM. If we take the cross section to be velocity dependent a — a Q (v/v a ) a , then the bound 
becomes [59"] 

°o/m x ( — ^ . ) <0.02cm 2 /g. (12) 



100 km/s 

It should be noted that, with a few exceptions outlined above, most simulations to date model DM self-interactions 
as a classical hard-sphere and thus assume the cross section is velocity independent, i.e. da/dfl — b 2 . By allowing the 
scattering cross-section to depend on inverse powers of v, the more restrictive bounds from clusters could in principle 
be avoided, while still allowing large cross sections at low velocities that could alter the structure of dwarf galaxies 
from the standard dark matter prediction. 

However, as the majority of the constraints are not derived from analyses that allowed for velocity-dependent cross 
sections, it is not completely clear whether we can apply these bounds to our work. One should bear in mind that 
the effects of large scattering are somewhat difficult to anticipate; for example, we expect larger cross sections to 
increase energy transfer rates, however if the mean free path becomes much smaller than the relevant length scales in 
the system then the medium becomes optically thick and transfer rates actually decrease. 

Since we cannot repeat a full iV-body simulation with the velocity dependent scattering cross section given by 
Eq. (|5j), we instead calculate the transfer cross section, averaged over the DM velocity distribution in the appropriate 
system, and compare to the corresponding bound on a/m x . The transfer cross section, rather than the unweighted 
cross section, is the parameter of interest since it measures the rate at which energy is transferred in the system. As 
outlined above, this is the quantity being constrained in the various limits, either directly or indirectly. In particular 
we calculate 

K) = J d 3 v, d 3 v 2 dCl - cos0)^ , (13) 

where in the galaxy's frame the WIMPs have velocity distribution 

Though we have, for completeness, listed a number of bounds in Table [I] we shall show that the one of most interest 
is the dwarf galaxy bound, as it comes from the system with the lowest dispersion velocity. For comparison, we shall 
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System 


w [km/ s] 


a/m x [cm 2 /g] 


References 


Bullet Cluster 


1000 


1.25 


snsa 


Galactic Evaporation 


1000 


0.3 


EH 


Elliptic Cluster 


1000 


0.02 


m 


Dwarf Evaporation 


100 


0.1* 


m 


Black Hole 


100 


0.02* 


m 


Mean Free Path 


44 - 2400 


0.01-0.6 


m 


Dwarf Galaxies 


10 


0.1 


m 



TABLE I: The systems we consider and the observational bound they place on DM self-scattering cross section. Entries marked 
with an asterisk * are velocity dependent bounds. For more details, see text. 



also consider in detail the bounds placed by elliptic clusters, as those represent the tightest bounds on a system with 
high velocity DM. 

In Fig. ^ we show (<7 tr )/m x as a function of m^, assuming a Maxwellian distribution with characteristic velocity 
t>o = 1000 km/s, which is approximately the value found in galaxy clusters. For values of greater than 0.5 GeV 
we include only I modes of zero and one, while for < 0.5 GeV, we include I < 5. As can be seen in Fig. [2j for our 
choice of a, rn^ and v, this is an acceptable trade-off between computational speed and accuracy. We also display the 
approximate solutions for the cross section and transfer cross section, as given by Eqs. ^ and (10), again integrating 
over a Maxwellian distribution for both incoming particles (the upper line is the approximate cross section, while the 
lower is a tI ). We can clearly see that for systems with velocity distributions centered around 1000 km/s no bounds 
on MeV-scale dark forces can be placed. 
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FIG. 5: (o"tr) /m x as a function of m$, assuming that m x — 500 GeV, and a = 0.01 (left) and a = 0.1 (right). A thermal velocity 



distribution Eq. (14 1 with dispersion vq — 1000 km/s = 3.3 x 10 c, characteristic of galaxy clusters, was used. Contributions 
from modes up to t = 5 are included in the exact numerical cross section for ra^ < 0.2 GeV, while only I < 1 are included 
above this mass. The approximate solutions from Eqs. (J9j and (10 1 are also shown (dashed red lines). 



However, dwarf galaxies, with velocity dispersions of ~ 10 km/s [61] . provide a non-trivial constraint. In Fig. [6] 
we show the velocity-averaged {atr)/m x as a function of m^, this time for a dwarf galaxy- appropriate value of 
v = 10 km/s. Again, both the exact numerical solution (with £ < 5 for all values of m^) and the approximate 
solutions are shown. Taking the upper bound on cr/m x to be the 0.1 cm 2 /g derived from dwarf galaxies, we can place 
a bound requiring 



m > 40 McV 



(15) 



for the larger value of a considered and slightly weaker (m^ > 30 MeV) for smaller a. Although clusters present a 
tighter bound on the scattering cross section, the characteristic velocity in these systems is far higher (Table [i]) and 
the stronger constraint comes from dwarf galaxies. 

A full simulation for the case of velocity dependent cross sections, as expected in models with a Sommerfeld 
enhancement, would improve on our estimate and we advocate strongly for it to be carried out. 
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FIG. 6: (ati) /m Y as a function of m,/,, assuming that m x = 500 GeV, and a = 0.01 (left) and a = 0.1 (right). A thermal velocity 



distribution Eq. (141 with dispersion Vo = 10 km/s = 3.3 x 10 _J c, characteristic of a dwarf galaxy, was used. Contributions 
from modes up to £ = 5 are included in the exact numerical cross section. The approximate solutions from Eqs. |9| and ( 10 1 
are also shown (red dashed lines). Regions of low m^> from the upper plots are shown in more detail in the lower plots. The 
limit from dwarf galaxy structure is 0.1 cm 2 /g. 



CONCLUSIONS 



An appealing explanation of the recent cosmic ray excesses is in terms of DM annihilation in our local galactic 
neighborhood, with the necessary increase in annihilation cross section (above that for a typical thermal WIMP) 
being due to a new dark force. The exchange of this dark force carrier </> - either a scalar or vector boson - leads to an 
attractive force between dark matter particles and consequently an increase at low velocities in their annihilation cross 
section: the Sommerfeld enhancement. To get a large enough enhancement the force must be sufficiently long-range 
and the force carrier is typically taken to be sub-GeV in mass. This is done for two reasons. The first is that the 
enhancement approximately scales as ~ am x /m^,, although it should be noted that for Yukawa potentials there are 
nearly massless bound states that lead to a far greater enhancement and a scan over parameter space shows that the 
density of these bound states is relatively high. The second is that the cosmic ray excesses are in leptonic channels 
and one way this might occur is if the DM annihilates to <fi particles which in turn decay to SM states. If the force 
carrier is sub-GeV in mass scale its decay products will be leptons. 

Solving the Schrodinger equation for the case of a Yukawa potential between the DM allows us to go beyond the 
Born approximation, resumming the effects of multiple light boson exchange, and to capture the effects of resonances 
in the potential. We have shown that the physics which causes the large enhancement in annihilation necessarily 
implies an increase in the self-scattering cross section. Requiring that this scattering enhancement does not cause 
significant energy transfer inside dark matter halos allows us to constrain the possible mass of the mediator ttia,, 
independent of the details of how it couples to the SM. There are various bounds on the self-scattering cross section 
coming from multiple systems. Approximately speaking, they amount to requiring less than one DM scattering per 
Hubble time, but by comparing iV-body simulations to observations this approximate bound can be improved upon 
in certain cases. 

Unfortunately many of these systems have very different DM velocity distributions and this hinders easy translation 
from a bound on one system to another, as most A^-body simulations assume a velocity independent hard sphere cross 
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section for DM scattering. We advocate strongly for TV-body simulations to be carried out which include the velocity 
dependence appropriate to the light mediator in the DM cross section. To allow comparison between these bounds and 
our velocity dependent scattering we calculate the transfer cross section averaged over the DM velocity distribution 
in two systems: a galaxy cluster and a dwarf galaxy. We find that galaxy clusters offer no interesting bound whereas 
dwarf galaxies require that > 40 MeV, with some sensitivity to the DM-mcdiator coupling. For mediators heavier 
than this bound there is still considerable Sommerfeld enhancement, both close to and away from resonances. This 
bound may well be improved upon by detailed simulations and or as further observations become available. 

Direct searches for a new U(l) gauge boson which would play the role of (f> have been suggested in e~e + collisions, 
and will be sensitive to of 0(1 GeV) [62] . Current bounds already exist for a range of and photon-0 mixing 
from T(3S) decay at BaBar [63) . measurements of e and fi anomalous magnetic moments [64] , the beam dump 
experiments at SLAC [5SJ [BS] and Fermilab 67], and dark photon searches at the Tevatron [BS]. However, these 
experiments, while probing the same mass range as Eq. (15) depend not only on m^, but also on its couplings to the 
SM. Our result is independent of the mixing with the photon, and in fact holds regardless of whether <j> is a U(l) 
gauge boson, a gauge boson of a non-abelian gauge group, or a scalar field. In that sense, the result is very robust. 

In conclusion, we have found that, for an experimentally interesting range of parameters, the proposed new 'dark 
forces' of interest in explaining results from both direct and indirect DM searches would have measurable impact on 
the structure of dwarf galaxy halos. As the enhancement of dark matter self-scattering does not depend on the details 
of the attractive force carrier beyond its mass and self-coupling, these bounds are applicable across a wide range 
of possible scenarios. Indeed, the increase in the scattering cross section could be present even in models where DM 
self-annihilation was forbidden by some symmetry. However, at this point, more work is necessary in order to be 
confident that the bounds on <r/m x , derived from Y-body simulation assuming a velocity independent cross section, 
can be applied in to situations where a tr depends on inverse powers of v. From the calculations provided in this paper, 
we believe that investigating this possibility is a worthwhile task, as it does seem that the effects are large enough to 
offer interesting constraints. 
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Note Added: During the completion of this work we became aware of a paper that addresses similar issues |38| . 
Though they consider a different astrophysical system we reach similar conclusions. 
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